arXiv: 1503.05194vl [astro-ph.GA] 17 Mar 2015 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 


Printed 19 March 2015 (MN MeX style file v2.2) 


Galaxy Pairs in the Sloan Digital Sky Survey - X: Does gas content 
alter star formation rate enhancement in galaxy interactions? 

Jillian M. Scudder^’^*, Sara L. Ellison^ Emmanuel Momjian^, Jessica L. Rosenberg‘S, 
Paul Torrey^’^’^, David R. Patton^, Derek Eertig"^, J. Trevor Mendel^ 

* Department of Physics and Astronomy, University of Victoria, Victoria, British Columbia, V8P lAl, Canada. 

^ Astronomy Centre, Department of Physics & Astronomy, University of Sussex, Brighton, BNl 9QH, England. 

^ NRAO, Domenici Science Operations Center, Socorro, NM, 87801, USA. 

^ School of Physics, Astronomy, and Computational Science, George Mason University, Fairfax, VA, 22030, USA. 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA. 

® Department of Physics, Kavli Institute for Astrophysics and Space Research, Massachusetts Imtitute of Technology, Cambridge, MA 02139, USA. 

^ TAPIR, Mailcode 350-17, California Institute of Technology, Pasadena, CA 91125, USA. 

* Department of Physics U Astronomy, Trent University, Peterborough, Ontario, K9J 7B8, Canada. 

^ Max-Planck-lnstitut fiir extraterrestrische Physik, 85748 Garching, Germany. 


19 March 2015 


ABSTRACT 

New spectral line observations, obtained with the Jansky Very Large Array (VLA), of a sam¬ 
ple of 34 galaxies in 17 close pairs are presented in this paper. The sample of galaxy pairs is 
selected to contain galaxies in close, major interactions (i.e., projected separations <30 
kpc, and mass ratios less extreme than 4;1), while still having a sufficiently large angular 
separation that the VLA can spatially resolve both galaxies in the pair. Of the 34 galaxies, 
17 are detected at > 3cr. We compare the Hi gas fraction of the galaxies with the triggered 
star formation present in that galaxy. When compared to the star formation rates (SFRs) of 
non-pair galaxies matched in mass, redshift, and local environment, we find that the star for¬ 
mation enhancement is weakly positively correlated (~ 2.5cr) with Hi gas fraction. In order 
to help understand the physical mechanisms driving this weak correlation, we also present 
results from a small suite of binary galaxy merger simulations with varying gas fractions. The 
simulated galaxies indicate that larger initial gas fractions are associated with lower levels of 
interaction-triggered star formation (relative to an identical galaxy in isolation), but also show 
that high gas fraction galaxies have higher absolute SFRs prior to an interaction. We show 
that when interaction-driven SFR enhancements are calculated relative to a galaxy with an 
average gas fraction for its stellar mass, the relationship between SFR and initial gas fraction 
dominates over the SFR enhancements driven by the interaction. Simulated galaxy interac¬ 
tions that are matched in stellar mass but not in gas fraction, like our VLA sample, yield the 
same general positive correlation between SFR enhancement and gas fraction that we observe. 

Key words: galaxies: interactions - radio lines: galaxies - galaxies: star formation - galaxies: 
starburst 


1 INTRODUCTION 


The interactions between galaxies are an ideal laboratory to study 
the response of galaxies to strong external perturbations. Galaxy- 
galaxy interactions very clearly drive strong changes in the per¬ 
turbed galaxies. Interacting galaxies are universally observed to 
have higher than average star formation rates (SFRs; e.g., [Li 


bas et al.|2003l|Alonso et al.|2004||Woods & Geller|2007l [Eliii 


el^a y20d8 2010 |Xu et al.|2010[ Darg et al.|2010 Stierwalt et al. 

201411, diluted nuclear gas-phase metallicities (e.g., Kewley et al. 
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2006 IMichel-Dansac et al.|2008||Ellison et al. 12008 [Scudder et al.| 


2012b^ , bluer than average nuclear colours (e.g., Patton et al.|201 1|(, 
and an overabundance of Active Galactic Nuclei (AGN; e.g., Koss 


et al.|201 lllSilverman et al.|201 11 

Ellison et al.|201 lllSabater et al. 

2013| Ellison et al.|2013 Satyapa 

et al.|2014| 

Khabiboulline et al. 


The combination of observational results with the analysis of 
simulations of galaxy interactions allows us to interpret these data 
as the many indirect signs of gas flowing to the centre of the galaxy 
in response to the interaction. These gas flows are understood theo¬ 
retically as the result of gravitationally driven tidal torques induced 
by the companion galaxy. These torques act upon misaligned stel- 
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lar and gaseous bars within the perturbed galaxies, produced as a 
result of the non-axisymmetric nature of the gravitational well, re¬ 
sulting from the presence of a companion galaxy. The misalignment 
between the gaseous and stellar hars allows the stellar component 
to effectively remove angular momentum from the gaseous com¬ 
ponent, the loss of which causes the gas to fall towards the centre 
of the galaxy, producing the flow needed to explain the observa- 


tions (e.g., |Mihos & Hemquist||1994| 

1996[|Barnes & Hernquist| 

19961 [Cox et al.|2006l|Di Matteo et al. 

2007l|Montuori et al. [20101 

Torrey et al.|201^. The gas inflow builds up a considerable reser- 


voir of gas near the centre of the galaxy, which can then be ef¬ 


ficiently converted into stars, fuelling a nuclear starburst (e.g., Di 


Matteo et al.|2007[[Montuori et al.|2010||Torrey et al.|2012[[M oreno| 

et al.|2014 1 . While galaxies in close pairs on average show higher 
than expected SFRs, the level of triggered star formation can vary 
widely from galaxy to galaxy, and from interaction to interaction 
( [Scudder et al.|[2012b^ . Many previous studies have attempted to 
isolate parameters which govern the strength of the star formation 
triggered within a galaxy. A few parameters are consistently found 
to modulate a galaxy’s response. 

The mass ratio between the two galaxies permits the strongest 
SFR enhancements when the two galaxies are approximately of 
equal mass. Unequal mass mergers are observed to be much less 
effective at triggering strong starbursts ( >^ods^Gellerj2007 El 


lison et al.|2008[|Xu et al.|2012||Lambas et al.|2012||Scudder et al. 


2012b| (, as the gravitational influence of a low mass companion 
should be weaker jCox et al.|2008||jret al.|2014) . The environment 
of the pair can also suppress triggered SFR if the galaxy is found 
within a particularly dense region (e.g., [Alonso et al.|2004|[2006) 
[Baldry et al.|2006[ [Ellison et al.|2010^ . Observational studies have 
shown that the strongest enhancements in SFR are visible at the 


smallest projected separations (r^; e.g., Lambas et al. 2003 Alonso 
[et al.[2004[rWoods & Geller|2007|[Ellison et al.|2008||scudder et ah 

[2012b[ [Patton et al.[[2013t . Unlike the simulations, observational 
results are limited to a snapshot view of any individual interaction; 
r,, serves as the best observational proxy for timescale and merger 
stage, properties which are predicted by theoretical studies to also 
modulate the SFR response of a galaxy interaction (e.g., |Torrey| 
[et al.[2012[ (. 

Theoretical works have also offered several additional param¬ 
eters which could alter the strength of a triggered starburst. For 
instance, the increasing prominence of a galactic bulge should sup¬ 
press triggered SFR (e.g., [Di Matteo et al.[2007[ (. Most simulations 
also agree that the orbital parameters imposed at the beginning of a 
simulation can have a significant effect upon the strength of the trig¬ 
gered star formation; the inclusion of a range of orbital parameters 
into a simulation suite introduces significant scatter, (e.g., Mihos[ 


et al.|1992[|Di Matteo et al.'|2Q07l[Cox et al.[2008[|Perez et al.|201 1 


Torrey et al.[2012[ l. Generally, any configuration or orbit of the two 
galaxies which allows for a particularly strong tidal force to act 
upon the two galaxies will result in a stronger triggered starburst. 
Flowever, the details of the orbital parameters of an interacting sys¬ 
tem are not readily constrained with observations, and therefore are 
generally left as a known source of scatter in the observations. 

The gas fraction of the interacting galaxies has also been sug¬ 
gested by the simulations to play a role in the triggering of SFR 
in galaxy pairs. Flowever, the impact of gas fraction has been less 
well-developed theoretically than the impact of mass ratio or orbital 
inclination. As such, theoretical works have not yet built a compre¬ 
hensive framework to understand the impact of varying gas fraction 
on the SFRs of galaxies in pairs. Different theoretical studies pro¬ 
vide windows into the physical processes that may be at play. Pre¬ 


vious works ( [Boumaud et al.|201 l[[Perez et al.|2011[ [Ferret et al.[ 
[2014[ ( have investigated the SFR response in high redshift interac¬ 
tions, constructed to have high gas fractions. Galaxies with high gas 
fractions are found to sustain strong, often turbulent, star formation 
well before the two simulated galaxies interact. The interpretation 
of these results varies; strong enhancement in the SFR is described 
by [Bournaud et al.'] ( [2011), Perez et al. ([2011^ found enhancement, 
but at a low level, and [Perret et al.[ ( [2014[ ( found no evidence for 
additional star formation due to the interaction. [Ferret et al.[j20l^ 
suggests that at sufficiently high SFRs, the galaxy may be able to 
‘saturate’ its star formation, unable to drive SFR to higher values 
even though the interaction causes a significant perturbation. 

Flowever, these high redshift analogue simulations do not ex¬ 
plicitly test the effect of varying the gas fraction within a galaxy 
on the SFR. This is undertaken by [Hopkins et al.[ ( [2009t , which 
suggests that the gas fraction of a galaxy is anti-correlated with its 
starburst efficiency. This anti-correlation arises due to the physics 
behind the origin of the gas flows in interacting galaxies. If the 
merger’s starburst strength is driven by the ability of the torque be¬ 
tween the stellar and gaseous bar to efficiently strip the gaseous 
bar of angular momentum, the relative strengths of those bars may 
overwhelm the importance of the volume of gas available. A high 
gas fraction galaxy would be very inefficient at driving gas flow, 
as a relatively weak stellar bar would exert a much smaller force 
upon a more massive gaseous bar. The weaker force would not be 
able to drive large quantities of gas to central regions of the galaxy, 
and the nuclear starburst which followed would be marginal. A low 
gas fraction galaxy would show the inverse effect; a very efficient 
torque of the small amount of gas within the galaxy, an effect also 
suggested by [Cox et al.[ ( [2008) . 

In spite of the predicted relationship between increasing gas 
fraction and decreasing triggered star formation within the galaxy 
interaction, [Di Matteo et al.[ ( [2007^ found that the orbital param¬ 
eters of an encounter drove far more scatter in the amount of star 
formation triggered in an interaction than variations in the gas frac¬ 
tion did. The large scatter found by [Di Matteo et al^p007[ ( implies 
that any correlation between SFR enhancement and gas fraction 
may be buried due to the impact of including varying orientations. 

1.1 The current work 

To date, there is no comprehensive study of H i gas in interacting 
galaxies which examines their individual SFR properties. There are 
a few studies which, while they contain no direct H i observations, 
interpret their results in the context of gas fraction. Studies of SFR 
triggering in high redshift galaxies (presumed to be at higher gas 
fractions) have found evidence for moderate SFR triggering jjogee 


et al. 2009[[^ baina et al.|2009[ |Wong et al.|201 l[[Whitaker et al. 

2014 ', but cf. Xu et al. ^2^^J7^which found no evidence for SFR 

enhancement. Historically, most direct studies of the gas content 
of interacting galaxies have been limited to detailed studies of one 
or two galaxies (e.g., [Smith! 199 1[[T9^ [Hibbard et al.|1994[[Du^ 


[et al.[1997[[Hibbard et al.[2()01[[Sengupta et al.[2013[l. 

There are several blind Hi surveys of the sky which could be 
used to investigate galaxy interactions, such as the H i Parkes All- 
Sky Survey (HIPASS; [Barnes et al.|200lf and the Arecibo Legacy 
Fast ALFA (Arecibo L-band Survey Array) survey (ALFALFA; 
[Giovanelli et al.|2005[ [Haynes et al.|2011[ l. ALFALFA has signif¬ 
icant overlap with the Sloan Digital Sky Survey (SDSS) footprint, 
and so could be very useful. A number of studies have made use 
of samples of galaxies cross matched between ALFALFA and the 
SDSS (e.g., [Catinella et al.[2012[[Stierwalt et al.[2014{[Fertig et al.[ 
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|2015[ |Ellison et al.|20l5) , though none of them have investigated 
the relationship between Alog(SFR) and gas fraction for individual 
galaxies. This is largely due to resolution limitations of single dish 
telescopes, which are unable to resolve individual galaxies in close 
pairs, where SFRs are most likely to be enhanced; gas fractions 
could therefore only be investigated as a pairwise property. 

In order to reach the spatial resolutions required to investi¬ 
gate the gas content of galaxies at small separations, interferomet¬ 
ric observations are required. We therefore take new observations 
with the NRAO’s Karl G. Jansky Very Farge Array (VFA) to ob¬ 
tain Hi masses for a sample of 34 galaxies in pairs. Our 17 galaxy 
pairs (34 individual galaxies) from the SDSS Data Release 7 (DR7; 
|Abazajian et al.|2009| l have pre-existing SFRs and stellar masses. 
With the pre-existing data products from the SDSS, we calculate 
gas fractions and determine whether a strong observational corre¬ 
lation exists between gas fraction and SFR enhancement. 

In Section|^ we present our sample selection criteria. In Sec¬ 
tion]^ we present the data reduction steps for the acquired H i data. 
In Section]^ we present the reduced data products from the VFA, 
along with the correlation between SFR enhancement and the gas 
content of the galaxies. In Section|^ we compare our observational 
results with a small suite of simulations. In Section]^ we discuss 
the results, and present our conclusions in Section]^ Throughout 
this work, we assume Q.m = 0.3, Qa = 0.7 and Ho = 70 km s“' 
Mpc“*. 


2 SAMPLE SELECTION 


Our master sample of galaxy pairs comes from the catalogue of 
|Patton et al.H20lf) . This catalogue contains > 23,000 galaxies in 
spectroscopically confirmed pairs from the SDSS Data Release 7 
(DR7) within Vp < 80 kpc, and with relative velocity differ¬ 
ences (Av) < 10,000 km s“'. In order to exclude galaxies which 
are likely to be projected pairs, previous works have traditionally 
imposed a limit on the Av of the pair somewhere between Av < 300 
km s“* (e.g., Scudder et al.||201^ i and Av < 500 km s“' (e.g., 
IFllison et al.| 2008^ . For the current work, we impose a limit of 
Av < 300 km s“* onto the galaxy pairs catalogue. 

We require that all potential target galaxies have well deter¬ 
mined stellar masses as determined by Mendel et al.|p0T4) and 
SFRs determined by |Brinchmann et ^ 2004^ for the SPSS DR7. 
Since the goal of this work is to examine the dependence of the 
interaction-triggered SFR enhancement on the gas fraction of the 
galaxy (which is defined as the H i mass relative to the stellar mass), 
these quantities are critical for our analysis. To ensure that the SFRs 
determined by [Brinchmann et af] ( |2004[ l are uncontaminated by 
any flux from an AGN, we also require thaf at least one galaxy 
in the pair be classified as star forming on the |Baldwin et al. J1^8H 
diagram, using the classification scheme described by |Kauffmann | 
|et al.| ( |2003| l. In order to be placed on the classification diagram, 
the S/N of the four required emission lines must be > l.cQ Only 
one galaxy in our sample of 17 pairs is classified as an AGN, and 
is paired with a galaxy which is classified as sfar forming. We note 
that all of our galaxies are classified as star forming by the slightly 
more permissive |Kewley et al.| ( |T001^ classification. The remainder 
of the sample are star forming galaxies. We note that this criterion 
requires that all galaxies be emission line galaxies. While this is not 


' Although the minimum S/N permitted is 1.0, all of the galaxies in the final 
sample have S/N > 3.0 in the four emission lines required for classification. 


Av 

Redshift 
Mass ratio 

Angular separation 


< 300 km s * 
0.01 < z < 0.06 

0.25< Ml!Ml <4.0 

< 30 kpc 
> 35 arcsec 


Table 1. Sample selection requirements for the galaxy pairs sample, in rel¬ 
ative velocity difference (Av), redshift (z), mass ratio, projected separation 
(Vp), and angular sky separation. 


an explicit cut on morphology, the majority of the sample will tend 
towards late-type morphologies. 

We also impose a lower limit on the redshift of our galaxy 
pairs of z > 0.01, as below that threshold, the deblending algorithm 
of the SDSS is prone to shredding galaxies into multiple objects. 
An upper redshift limit of z < 0.06 was applied in order to prevent 
the exposure times for the VLA from becoming prohibitively large. 

In order to assemble a sample of galaxy pairs which is most 
likely to be showing signs of triggered star formation, we use the 
results of |Scudder^eWlJj 201^M to impose additional criteria on 
our sample. Scudder et al. {lQ\2h\ found that galaxies in pairs are 
most likely to be strongly enhanced when found at small Vp from 
their companion (i.e., rp < 30 kpc), and in nearly equal mass 
mergers (0.33 < MhostlMcompanion < 3.0). We accordingly impose 
the criterion that the galaxy pairs must be at separations < 30 
kpc. In order to prevent the sample from becoming unnecessarily 
small, we opt to relax the definition of a ‘major merger’ slightly and 
impose the criterion 0.25 < Mi,osiIMcompanion < 4.(j^ This slightly 
more permissive definition has been used in the literature in the past 
as an alternate criterion to select nearly equal mass mergers (e.g., 
[Woods & Geller|2007[|Lotz et al.|2010t . 

Our final criterion was to eliminate the galaxy pairs which 
were too closely positioned on the sky for the VLA’s C config¬ 
uration to be able to spatially resolve the two galaxies. The syn¬ 
thesised beam width of C configuration is ~ 14 arcsec, assuming 
uniform weighting. We therefore required that the central locations 
of the two galaxies be separated by > 35 arcsec (or 2.5 times the 
synthesized beam width). In combination with the criterion that the 
galaxies be at physical separations < 30 kpc, this biases our 
final sample towards the low redshift end of the permitted range. 

The full set of sample selection requirements is summarised 
in Table [T] and results in a sample of 17 galaxy pairs, containing 
34 individual galaxies. The galaxies in our sample are presented in 
Table 1^ which lists the SDSS Object ID (objid), the between 
pairs in kpc, the velocity difference between the pair in km s“*, 
the mass ratio between the two galaxies, and the stellar masses of 
each galaxy. 


3 DATA ACQUISITION 

Of the 17 galaxy pairs in our final sample, we were awarded VLA C 
configuration time for 12 galaxy pairs (24 galaxies) under Proposal 


^ When this sample was defined, the stellar masses (and by extension, the 
mass ratios) were the same as those used in jScudder et al.|j2012b|. Since 
the acquisition of the data, the mass estimates have been refined by|Mendel| 
jetal.jpOTT) . In general, this does not appreciably change the stellar masses, 
but for two pairs in our sample, the stellar masses were changed by a suffi¬ 
cient amount to push their mass ratio past the 4:1 limit. They have not been 
excluded from the sample. 
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SDSS ObjID 


Av 

Mass ratio 

log M, 



(kms-‘) 


(Mo) 

587726033341776175 

21.51 

36 

3.66 

9.65 

587726033341776191 




10.21 

587739609695453284 

22.27 

65 

2.50 

9.77 

587739609695453281 




10.17 

587744873717563559 

20.18 

122 

4.26 

9.14 

587744873717563471 




8.51 

587727179536859247 

15.45 

111 

3.70 

8.90 

587727179536859227 




8.33 

587729160043757697 

22.03 

48 

2.90 

8.61 

587729160043757707 




9.08 

587741489815027774 

26.78 

91 

1.12 

9.83 

587741489815028146 




9.78 

587739303684866183 

20.2 

293 

3.03 

8.61 

587739303684866173 




9.09 

587726033308680234 

28.57 

71 

1.27 

9.65 

587726033308680320 




9.75 

588018056204780081 

29.46 

208 

4.64 

9.81 

588018056204780049 




10.48 

587729158970867777 

16.77 

47 

1.31 

9.87 

587729158970867792 




9.76 

587733605328093368 

27.7 

130 

1.15 

9.18 

587733605328093256 




9.12 

588848899908370674 

23.46 

16 

2.20 

9.19 

588848899908370505 




9.53 

588017605758025795 

15.2 

49 

1.09 

8.46 

588017605758025732 




8.50 

588017702411763744 

22.04 

68 

3.66 

9.64 

588017702411763872 




9.08 

587727178473930875 

27.87 

73 

1.24 

8.94 

587727178473930886 




8.85 

587742901789589569 

19.55 

118 

3.59 

9.73 

587742901789589575 




9.17 

588023670245949622 

19.39 

53 

1.63 

8.35 

588023670245949625 




8.56 


Table 2. SDSS objids, projected separations, velocity differences, mass ra¬ 
tios, and stellar masses for each of the 17 galaxy pairs in the final sample. 

ID 12A-061, with a mixture of Priority A, B, and C time. Director’s 
Discretionary Time was successfully obtained during Cycle 13 A to 
obtain VLA C configuration data for the remaining 5 galaxy pairs 
(Proposal 13A-537). Of the awarded time, almost all of the schedul¬ 
ing blocks were observed; only one target (588023670245949622 
& 588023670245949625) was missing on-source time when the 
telescope ceased C configuration observations. In this section, we 
briefly summarise the setup of these observations. 

3.1 Telescope setup 

In both proposals 12A-061 and 13A-537, we used the C config¬ 
uration of the VLA, which has a maximum baseline between an- 


Receiver: L Band 1000 - 2000 MHz 

Rest Frequency 1420.40575 MHz 

Sub-band Bandwidth 4.0 MHz 

No. of poln. products: 2.0 

No. of channels/poln product 128 

Channel Width (Frequency) 15.625 kHz 

Channel Width (Velocity) 3.30kms“*@ 1420 MHz 


Table 3. Project 12A-061 telescope setup 

tennas of 3.4 km. C configuration is one of the more compact of 
the VLA’s configurations, allowing it to be more sensitive to large 
scale structure on the sky, relative to the more extended configu¬ 
rations of the array (A & B). However, this large-scale sensitivity 
comes at the expense of angular resolution. C configuration has a 
synthesised beamwidth of 14 arcsec with a sensitivity to large scale 
features up to 970 arcsec in size. At the typical redshift of our sam¬ 
ple (z ~ 0.02), this large angular scale sensitivity corresponds to a 
physical size of ~ 390 kpc. The large synthesised beam of the 
array in this configuration means that in general, the Hi flux of a 
galaxy will be contained within a single beam, seen roughly as a 
point source by the array. As the galaxies are well separated both 
spatially and in velocity, there should not be contamination from 
the companion galaxy. 

For the data taken in 12A-061, we used the dual polarisation 
mode of the WIDAR correlator in Open Shared Risk Observing 
(OSRO) program. The data were taken between February 05, 2012, 
and March 18, 2012. The full spectral coverage was divided into 
2 sets of 8 subbands, where each subband in one set was offset 
from the second set by half the width of the subband. This resulted 
in high signal to noise coverage across the entire frequency range 
probed. The correlator integration time was 5 seconds. A summary 
of the telescope setup is presented in Table[^ 

The data acquired as part of project 13A-537 are also in dual 
polarisation mode, but the Open Shared Risk Observing mode used 
in 12A-061 had been retired as the telescope exited commissioning. 
The data acquired as a part of 13A-537 were observed between July 
26, 2013 and September 02, 2013. Instead of using multiple sub¬ 
bands offset from each other by half the width of the subband, in 
13A-537, a single 16 MHz wide subband was used, covering the 
entire frequency range. The 13A-537 observations were carried out 
at the same frequency resolution and time resolution as the data 
from 12A-061. A summary of the data acquisition is listed in Table 
this Table contains the SDSS objids, the redshifts of each galaxy, 
the r-band bulge to total ratio of the galaxy ( [Simard et al.|20lT| ), 
the coordinates in Right Ascension and Declination, the proposal 
ID under which the observations were taken, the date of the obser¬ 
vation, and the duration of the observations. 


3.2 Data reduction 

The VLA data were reduced and calibrated using the Common As¬ 
tronomy Software Applications (CASA ; |McMullin et al.|2007) and 
the Astronomical Image Processing System (AIPS; |Greisen l2003l > 
packages of the National Radio Astronomy Observatory (NRAO). 
The imaging and deconvolution were performed in CASA for all 
data sets. 

Following standard procedures, the data were flagged for in¬ 
strumental issues, radio frequency interference, and nonfunctional 
antennas, along with any other major problems noted in the observ- 
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ing logs. Bandpass calibration, complex gain calibration, and flux 
scaling were performed on each data set. 

The data were converted from the Fourier domain into the im¬ 
age plane using CASA’s clean task, and the line-free channels in 
the data cube were used to model and subtract off the continuum 
emission with uvcontsub. clean was then used to create image 
cubes. For those data sets in Table with multiple observations, 
CLEAN was used at this stage to combine the UV (Fourier) data into 
a single image cube. All images were created at a resolution of 4 
arc sec to a pixel. 

For each data cube, a spectrum was extracted for each galaxy 
in the pair. The optimal spectral region was selected, using the RA 
and Declination from the SDSS as a starting point. If necessary, the 
image cube was smoothed to coarser velocity resolution. Once the 
optimal spatial box and velocity resolution were determined, the 
spectrum for that galaxy was extracted in units of flux density per 
optical km s“*. 


4 DATA PRODUCTS & ANALYSIS 

In this section we describe the data products used for the analysis 
in the remainder of this paper. This includes the calculation both of 
the HI masses and of the upper limits to the H i masses, along with 
a more detailed description of the data products from the SDSS for 
our sample of galaxies in pairs. 


4.1 Hi masses 

For each galaxy in our sample. Hi masses were calculated using 
the standard equation: 

Mhi = (2.36 X 10^) J S(v) dv (1) 

S (v) is the emission flux in Jy km s“*. The complete integral 
f S (v) dv represents the integrated line flux over the width of the 
line. D is the distance to the object in Mpc. For this work, we calcu¬ 
late D as D = czjHo, which is applicable for low redshift objects, 
such as those contained within our sample. 

The significance of each H i mass was determined through a 
peak/RMS calculation. The peak value of the emission line is di¬ 
vided by the RMS of the spectrum outside of the velocity range of 
the galaxy’s emission. In Figure[T] we show two of the 17 acquired 
spectra. In the right hand side, we show the spectra for both galaxies 
in the pair, with the blue shaded region being the range of the spec¬ 
trum used for peak/RMS calculations. The vertical dashed line indi¬ 
cates the systemic velocity of the galaxy, with the grayscale shading 
indicating the uncertainty on the SDSS velocit)|^(±30 km s“'). The 
horizontal dotted line indicates the zero line. On the left hand side 
of Figure[T]we show the SDSS thumbnail image of the galaxy pair. 
The remaining 16 galaxy pairs in the sample are shown in the Ap¬ 
pendix (Figures [A2| through | A1 7| l. All calculated Hi masses, with 
their peak/RMS signal to noise estimates, are found in Table[^ 

Of the 34 galaxies in the sample, 17 are detected at S/N > 3.0. 
For galaxies whose spectrum has a marginal detection (peak/RMS 
S/N < 3.0), an upper limit is determined using the bounds of the 


^ The uncertainty on the SDSS redshift is an average value for the Main 
Galaxy Sample, as quoted at the following URL: http: //classic. sdss. 
org/dr7/ 


marginal spectral feature. If no marginal detection is found, a ve¬ 
locity width of Av= 300 km s“*, centred around the systemic ve¬ 
locity of the galaxy, is adopted. We use upper limit values for the 
remainder of the analysis for all galaxies which were detected at 
S/N < 3. Upper limit calculations use the same velocity width used 
for the non-detection, and assumes 3x the RMS value to be present 
across the entire velocity width. This integrated flux value is then 
input into Equation [T] and the output is taken as a S/N = 3 upper 
limit. These upper limits are included in Table 


4.1.1 Gas Fractions 


Our estimation of the gas fraction uses the following expression, 
following previous observational works: 


f - 

M, 


( 2 ) 


This is distinct from the calculation of most theoretical works, 
which prefer the form of: 


/sos — 




M, + 


(3) 


The calculated gas fractions are presented in Table We cal¬ 
culate gas fractions for all galaxies in the sample, including the 
upper limits, using the observational definition described in Equa¬ 
tion!^ For consistency with the observational results presented in 
this paper, we have presented all theoretical values using the form 
which is more directly comparable with the observational 
values, as the observations do not have total gas masses available. 

In order to calculate the uncertainties on any individual H i de¬ 
tection in our sample, for those galaxies with peak/RMS S/N > 3.0, 
we calculate a Icr error in a similar way to the calculation of the up¬ 
per limits. Instead of 3x the RMS value of the spectrum outside the 
spectral line region, we take lx the RMS outside the emission re¬ 
gion and calculate the equivalent H i mass associated with that flux 
density over the velocity width of the detected emission feature. 


4.2 SFR enhancements 

As our sample of galaxy pairs was selected to be part of the SDSS 
DR7, each galaxy in the sample has a large amount of ancillary 
data available, which includes stellar masses ( [Mendel et al.|201^ 
calculated from the updated photometry of jSimard et al. [201 1[ ( and 
SFRs ( [Brinchmann et al.|2004| l. In order to correlate the gas frac¬ 
tions with the star formation enhancement in each galaxy, we now 
must define a metric which quantifies the amount of star formation 
triggered in each galaxy due to their interaction. 


4.2.1 SDSS Control Sample 


In order to estimate the change to the SFR of a galaxy in the pairs 
sample (described in Sectionrelative to what is ‘expected’, we 
construct a robust control sample of galaxies which are not found 
in pairs. These control galaxies are selected to be tightly matched 
to the pairs sample in stellar mass, redshift, and local density. Local 
density is defined as in previous works (e.g., [Ellison et al.|2013|l as: 


d„ is the distance to the nearest neighbour within a redshift 
slice of ±1000 km s“’ around the redshift of the target galaxy; for 
this work, ii = 5. E5 is then normalised to the median value within 
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SDSS ObjID 

Redshift (e) 

B/T, 

RA (J2000) 

Dec (J2000) 

Proposal ID 

Date Observed 

Duration 

587726033341776175 

0.027775 

0.0 

14:58:21.64 

±02:32:48.94 

12A-061 

Feb 06 2012 

30 min 

587726033341776191 

0.027899 

0.0 

14:58:23.24 

±02:32:18.82 




587739609695453284 

0.027482 

0.0 

11:18:11.21 

±31:59:16.32 

12A-061 

Feb 22 2012 

30 min 

587739609695453281 

0.02726 

0.238 

11:18:10.97 

±31:58:35.93 




587744873717563559 

0.019347 

0.0 

09:01:30.48 

±12:39:43.63 

12A-061 

Feb 20 2012 

1 h 

587744873717563471 

0.019761 

0.223 

09:01:33.86 

±12:39:31.67 




587727179536859247 

0.018063 

0.0 

02:06:21.33 

-08:52:19.26 

12A-061 

Mar 17 2012 

1 h 30 min 

587727179536859227 

0.018439 

0.071 

02:06:20.21 

-08:52:57.50 




587729160043757697 

0.023771 

0.997 

13:32:02.01 

±05:32:12.73 

12A-061 

Feb 07 2012 

30 min 

587729160043757707 

0.023608 

0.0 

13:32:04.62 

±05:32:37.14 




587741489815027774 

0.03761 

0.0 

08:05:53.02 

±14:00:06.34 

12A-061 

Feb 06 2012 

30 min 

587741489815028146 

0.037924 

0.039 

08:05:55.42 

±13:59:59.04 




587739303684866183 

0.017681 

0.417 

13:12:07.51 

±34:16:13.78 

12A-061 

Feb 07 2012 

30 min 

587739303684866173 

0.016686 

0.18 

13:12:05.26 

±34:15:23.16 




587726033308680234 

0.031671 

0.752 

09:55:22.79 

±02:34:14.25 

12A-061 

Feb 09 2012 

30 min 

587726033308680320 

0.031426 

0.0 

09:55:23.34 

±02:34:58.81 




588018056204780081 

0.03083 

0.047 

17:05:01.61 

±23:09:26.43 

12A-061 

Feb 05 2012 

30 min 

588018056204780049 

0.030114 

0.072 

17:04:59.9 

±23:10:08.57 




587729158970867777 

0.023125 

0.0 

13:39:27.31 

±04:30:38.03 

12A-061 

Feb 07 2012 

30 min 

587729158970867792 

0.022964 

0.495 

13:39:29.71 

±04:30:40.95 




587733605328093368 

0.026354 

0.045 

13:59:27.75 

±58:19:33.72 

12A-061 

Feb 06 2012 

30 min 

587733605328093256 

0.025908 

0.192 

13:59:31.93 

±58:18:52.55 




588848899908370674 

0.025371 

0.0 

11:20:57.36 

-00:05:49.48 

12A-061 

Feb 10 2012 

30 min 

588848899908370505 

0.025317 

0.0 

11:20:57.06 

-00:05:03.73 




588017605758025795 

0.019123 

0.063 

12:18:38.72 

±45:46:28.13 

13A-537 

Aug 62013 

30 min 

588017605758025732 

0.019291 

0.224 

12:18:39.14 

±45:47:06.90 


Aug 09 2013 

30 min 







Aug 10 2013 

30 min 

588017702411763744 

0.029248 

0.154 

14:48:19.69 

±09:07:02.14 

13A-537 

Jul 26 2013 

30 min 

588017702411763872 

0.029015 

0.022 

14:48:17.59 

±09:07:23.48 




587727178473930875 

0.033197 

0.411 

03:45:43.95 

-06:54:46.02 

13A-537 

Aug 18 2013 

30 min 

587727178473930886 

0.032945 

0.205 

03:45:46.7 

-06:54:56.24 


Aug 19 2013 

30 min 







Aug 27 2013 

30 min 

587742901789589569 

0.025918 

0.906 

12:39:35.85 

±16:35:16.14 

13A-537 

Aug 13 2013 

30 min 

587742901789589575 

0.026322 

0.0 

12:39:38.35 

±16:35:06.34 




588023670245949622 

0.021097 

0.0 

12:20:37.91 

±20:40:21.19 

13A-537 

Aug 10 2013 

30 min 

588023670245949625 

0.021278 

0.655 

12:20:39.01 

±20:39:38.69 


Aug 13 2013 

30 min 







Sep 02 2013 

30 min 


Table 4. SDSS objids, redshifts, r-band bulge to total ratio, RA & Declination, proposal ID, date observed, and the duration of the observation block, for all 
galaxy pairs in the sample (34 galaxies in 17 pairs). 


a redshift range of ±0.01 to produce the normalised density (dj). 
is used to constrain the local density of the pairs and their control 
galaxies. 

Matching a control sample in stellar mass, redshift, and lo¬ 
cal density has been found to eliminate the majority of sources of 
bias between the two samples ( [Perez et al.|2009^ . Any remaining 
systematic issues within the data should be equally present in the 
pair sample and the control sample; any measurements of quantities 
which are relative between the two samples should not be affected. 

To be considered in the matching process, each potential con¬ 
trol galaxy must be classified as star forming on the BPT diagram 


according to the [Kauffmann et al.j p003| l classification scheme. 
This criterion, as for the pairs sample, requires that the S/N of the 
necessary emission lines be > 1.0, and excludes AGN and non¬ 
emission line galaxies. The control galaxy must also not be in¬ 
cluded in the master pairs catalogue of jPatton et al.| ( [^l 1[ (, which 
removes all galaxies with a spectroscopic companion at rp< 80 
/jyiJ kpc, and within Av < 10,000 km s“'. To further limit the spu¬ 
rious inclusion of galaxy pairs in the control sample which are 
missed due to spectroscopic incompleteness, we also require that 
the galaxy have a merger vote fraction = 0, as determined by 
the Galaxy Zoo participants ( jParg et al.|2010[ (. [Darg et ^ [2010[ ( 
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Figure 1. Left panel: SDSS thumbnail of galaxies 588018056204780081 (lower left) & 588018056204780049 (upper right). The image is 250 arcsec to a 
side. Right panels: Flux density vs. optical velocity for galaxy pair 588018056204780081 (upper) & 588018056204780049 (lower). The vertical dashed line 
indicates the systemic velocity of the galaxy, with the gray shaded region indicating the uncertainty on the SDSS velocity (±30 km s“*) The horizontal dotted 
line indicates the line of zero flux. The blue shaded region indicates the region of the spectra used for signal to noise and gas mass calculations. The upper 
spectrum has a peak/RMS S/N of 3.19, while the lower has a S/N of 1.07. The SDSS thumbnails and associated spectra for the remaining 16 galaxy pairs can 
be found in the Appendix. 


use a cut in the vote fraction of /,„ = 0.4 (i.e., 40% of Galaxy Zoo 
participants who inspected that particular galaxy classified it as an 
interacting system) to define a pairs sample. In order to be conser¬ 
vative in our exclusion of galaxy pairs, we use the limit of /„, = 0 
to ensure that interacting systems which were contentious or mis- 
classified by the Galaxy Zoo participants are also excluded. 

From this sample of potential control galaxies, we match each 
of the 34 galaxies in the VLA sample to a set of control galaxies 
which meet all of the criteria outlined above. Following the work 
of |Ellison et al.| ( f2013^ , control galaxies are matched to a paired 
galaxy if they fall within a redshift tolerance of Az=0.005, a stellar 
mass tolerance Alog M,=0.1 dex, and a local density tolerance of 
Ads = 0.1, relative to that of the galaxy in the pair. Typically, each 
pair galaxy is matched to > 50 control galaxies. 

Previous works (e.g., |Scudder et al.|[^12a|b[ |Ellison et al.| 
|2013| l have defined a quantity dubbed an ‘offset’, denoted as 
Alog(SFR). This offset value quantifies changes in SFR relative to 
the SFR - stellar mass relation of galaxies. We follow the proce¬ 
dure of |Ellison et al.| ( [2013| l in defining the offset value to be the 
difference befween the SFR of the pair galaxy and the median SFR 
of the non-pair control galaxies. 

Alog(SFR) = log (SFRp,„) - (log (SFR„„,„,,)) (5) 

Positive values of Alog(SFR) indicate that the paired galaxy 
has a higher SFR than its matched controls, while negative values 
indicate that the galaxy has a lower SFR than its controls. We cal¬ 
culate offsets for the total, aperture corrected SFRs (|Brinchmann| 
|et al.|2004{ [Salim et al.|2007| l, which provide the enhancement in 
the SFR across the entire disk. 

We obtain Alog(SFR) values for each of the galaxies in our 
pairs sample, which can then be individually plotted against the gas 
fractions calculated in Section [4.1.1| The median uncertainty in the 
SFR values among all star forming galaxies in the SDSS is 0.09 dex 


( jBrinchmann et al.|2004| ); we adopt this value as the uncertainty in 
the Alog(SFR) across the pairs sample. 


4.3 Correlating H i mass with Alog(SFR) 

Figure 1^ shows the relation between gas fraction and Alog(SFR), 
for those galaxies in the sample which were detected at S/N > 3.0. 
The colours and greek letter of each data point matches the two 
galaxies in a pair, with each pair plotted as a unique colour. Error 
estimates in both gas fraction and Alog(SFR) are shown in grey. 
The full sample is presented in Figure [AT| Galaxies with H i S/N < 
3.0 (non-detections) are plotted with their 3cr upper limits, denoted 
by smaller triangles. In both figures, it appears that there is a rough 
trend of the highest gas fractions being found in galaxies which are 
most strongly starbursting. 

To test the significance of this relationship between the total 
Alog(SFR) and the H i gas fraction, we make use of the Spearman 
Rank correlation test. This correlation statistic produces a coeffi¬ 
cient of 1.0 if the gas fraction of a galaxy always increases as it 
progresses to higher Alog(SFR), regardless of the shape of this cor¬ 
relation. Other correlation metrics typically test for linear correla¬ 
tions, which may be less applicable in this instance. A Spearman 
Rank correlation test of Figure|^produces a correlation coefficient 
of 0.5637. The Spearman Rank correlation test also outputs a null 
probability, which indicates the probability of no correlation, given 
the sample size. The probability of the data being uncorrelated pro¬ 
duced by this test is 0.0184. This excludes the null hypothesis at 
slightly less than 2.5(T. 

In addition to the Spearman Rank correlation test, we also run 
a least-squares fitting routine on the data in order to acquire the an 
estimate of the slope of the data in order to quantify the correlation. 
We use a line of the following form: 

log(M/;,/M,) = A*Alog(SFR)+S (6) 
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8 Scudder et al. 


SDSS ObjID 

Alog(SFR) 

log M. (Mo) 

log Mhi (Mo) 

(Mm/M.) 

Sf^(pk/RMS) 

587726033341776175 

0.25 

9.65 

<9.32 

<0.47 

2.85 

587726033341776191 

-0.11 

10.21 

8.68 

0.03 

3.43 

587739609695453284 

0.24 

9.77 

9.19 

0.26 

4.08 

587739609695453281 

-0.42 

10.17 

9.10 

0.09 

3.69 

587744873717563559 

0.001 

9.14 

8.48 

0.22 

5.35 

587744873717563471 

0.16 

8.51 

8.32 

0.64 

3.51 

587727179536859247 

-0.05 

8.90 

9.48 

3.82 

11.40 

587727179536859227 

0.02 

8.33 

< 8.78 

<2.81 

2.02 

587729160043757697 

0.60 

8.61 

8.64 

1.07 

3.77 

587729160043757707 

0.51 

9.08 

8.70 

0.42 

3.16 

587741489815027774 

0.27 

9.83 

<9.92 

< 1.23 

2.70 

587741489815028146 

0.54 

9.78 

<9.84 

< 1.15 

2.23 

587739303684866183 

-0.18 

8.61 

<9.19 

<3.77 

2.35 

587739303684866173 

-0.12 

9.09 

8.73 

0.43 

3.57 

587726033308680234 

0.32 

9.65 

<9.94 

< 1.97 

1.97 

587726033308680320 

0.41 

9.75 

<9.76 

< 1.02 

2.00 

588018056204780081 

0.18 

9.81 

9.32 

0.32 

3.19 

588018056204780049 

0.34 

10.48 

< 8.64 

<0.01 

2.17 

587729158970867777 

0.02 

9.87 

9.71 

0.69 

4.09 

587729158970867792 

-1.10 

9.76 

9.12 

0.23 

3.70 

587733605328093368 

-0.16 

9.18 

< 8.84 

<0.46 

2.26 

587733605328093256 

0.28 

9.12 

9.66 

3.53 

4.04 

588848899908370674 

0.22 

9.19 

9.60 

2.57 

4.54 

588848899908370505 

0.15 

9.53 

8.57 

0.11 

3.26 

588017605758025795 

-0.26 

8.46 

<9.17 

<5.64 

1.93 

588017605758025732 

0.70 

8.50 

9.22 

5.28 

3.62 

588017702411763744 

0.68 

9.64 

< 10.25 

<4.07 

2.67 

588017702411763872 

0.082 

9.08 

<9.96 

<7.68 

2.79 

587727178473930875 

0.67 

8.94 

< 10.02 

< 11.90 

2.62 

587727178473930886 

0.14 

8.85 

<9.50 

<4.46 

2.24 

587742901789589569 

0.54 

9.73 

9.54 

0.66 

3.20 

587742901789589575 

0.16 

9.17 

<9.99 

<6.53 

2.70 

588023670245949622 

0.16 

8.35 

<9.01 

<4.51 

1.61 

588023670245949625 

-0.31 

8.56 

<9.39 

<6.68 

2.73 


Table 5. SDSS ObjID, Alog(SFR), stellar mass, Hi gas mass, gas fraction, and the peak/RMS S/N for each of the galaxies in our sample. 3ct upper limits are 
included for all galaxies which were detected at S/N < 3. 


The fitting procedure incorporates the errors in both gas frac¬ 
tion and in Alog(SFR) as weights on each point during the fit¬ 
ting. The resulting line of best fit coefficients are A= 1.02 ± 0.42, 
B=-0.16 ± 0.17, where the uncertainties are the standard errors 
calculated for each coefficient. This fit is plotted in Figure]^ as the 
dashed grey line. Flowever, we find that the yellow outlying poinQ 
at Alog(SFR)= -1.0 skews the fit towards a shallower slope and 
does not produce a representative fit to the majority of the points. In 


Even though this galaxy has a very low overall Alog(SFR), it is still a 
star forming galaxy with significant (S/N > 3.0) emission lines. This galaxy 
has low SFR in the region outside the SDSS fibre, so even though it has 
significant SFR within the fibre region, the overall SFR is shifted to a low 
value. 


order to more accurately fit the remaining points, we exclude this 
data point from the fitting procedure. (The exclusion of this data 
point still results in a Spearman rank null correlation exclusion of 
> 2cr.) The best fit relation to the data if the outlier is excluded has 
coefficients A = 2.39 ± 0.82 and B = -0.59 ± 0.26. This relation is 
plotted in Figure|^as the dark grey line. 

To test the sensitivity of the Spearman rank correlation test to 
the distribution of our data points given the small number statistics, 
we also perform a jackknife resampling procedure. We remove one 
galaxy from the sample at a time, and re-perform the Spearman 
rank test on the 16 remaining galaxies, recording the null probabil¬ 
ity for each subsample. This is repeated until each galaxy has been 
removed from the sample once. We find that the probability of no 
correlation is remarkably robust to the removal of a single point; 
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Total A log(SFR) 

Figure 2. The log of the gas fraction vs. Alog(SFR) for galaxies with 
S/N > 3.0. Galaxies are matched in colour and hy greek letter if they 
are in the same pair. The dashed grey line indicates the results of a least- 
squares fit to the observations, using the horizontal and vertical error bars 
as weights on the least squares fit, and is described by logiMni/M,) = 
1.02xAlog(SFR)-0.16. The solid dark grey line shows the results of the 
same fitting routine, but excludes the point plotted at Alog(SFR) = -1.0 
(labeled with k), and is described by log(M/r//M,)=2.39xAlog(SFR)-0.59. 


the median null probability of the subsamples is 0.0235, which is 
still well above The probability of the data being uncorrelated, 
as calculated by the Spearman rank test, is therefore unlikely to be 
driven by single outliers. 


5 SIMULATIONS 

To develop a more complete physical picture of what could drive 
the observed ~ 2.5cr correlation between increasing gas fraction 
and increasing Alog(SFR), we turn to simulations. The existing 
theoretical work is not easily comparable to observations, as gas 
fractions are often measured in unobservable units, or compared at 
different time steps. As a comparison to theoretical work may illu¬ 
minate the physical mechanisms behind the observational results, 
we develop a small suite of simulations with the goal of comparing 
to observations. 


5.1 Simulation suite design 

Our suite is based on [Torrey et al.| ( |2012| ), and a complete de¬ 
scription of the model setup and parameters can be found in that 
work. Briefly, the simulations of |Torrey et al.| ( |2()T^ use GADGET-2 
( |Springel|2005| ) to construct a set of N-body/smooth particle hydro- 
dynamic models of binary galaxy mergers. These mergers include 
radiative gas cooling ( Katz^e£^J]^96^^^£to formation and associ¬ 
ated feedback ( [Springel & Hemquist|2003) , and chemical enrich¬ 
ment. The SFR of the pair can be tracked either across the entire 
disk or within a smaller region at the centre of the disk; to compare 
to the current observations, we use the SFR across the entire disk. 


The galaxies in the simulations are stable and do not develop a bar 
for at least 2 Gyr if they are evolved in isolation. 

In order to compare the results found here with the simula¬ 
tions, we construct a suite of binary mergers which vary only in 
gas fraction. From |Torrey et al.]j2012^ , we use orbit ‘e’. The sim¬ 
ulated galaxy models have initial stellar masses between 10.26 < 
log Mq < 10.52. For each gas fraction, the galaxy model is merged 
with an identical copy of itself. We then track log(SFR) along with 
gas fraction for the duration of the simulation. The suite uses initial 
gas fractions of = 0.04, 0.09, 0.19, 0.32, 0.47, 0.67, 0.92, 

1.27, and 1.78, for a total of 9 runs. 


5.2 Results from the simulatiou suite 

In the top panel of Figure we show the time evolution of the 
physical (3D) separation between the two galaxies. Each colour in¬ 
dicates a unique initial gas fraction, labeled in the legend, progress¬ 
ing from navy for the most gas poor galaxy interaction {MgasIM, = 
0.04) to pale blue for the most gas rich - 1.78). All of 

the galaxies have the same initial orbit; the only deviations are near 
coalescence, where slight discrepancies arise due to differences in 
dynamical friction. 

In the second panel, we show the log of the average SFR of the 
2 interacting galaxies, as a function of time. For equal mass galax¬ 
ies, the SFR response to the interaction should be approximately 
equal (e.g., |Torrey et al.|2012| ). This panel indicates that higher gas 
fraction galaxies have systematically higher SFRs across the entire 
interaction (and prior to any encounter) than the low gas fraction 
galaxies. 

In the third panel, we show the total Alog(SFR) for the galaxy 
in the interaction. Alog(SFR) is calculated in a very similar man¬ 
ner to the observations. For each time step in the simulatiorj^ 
we take the log(SFR) of the interacting galaxy, and subtract the 
log(SFR) of the same galaxy model in isolation at that same time 
step. Prior to their first encounter, the two galaxies have the same 
SFR (Alog(SFR) = 0.0). The lowest gas fraction galaxies show 
rapid strong enhancement (relative to their isolated counterpart) 
after a close encounter, with the Alog(SFR) remaining high and 
nearly constant for the duration of the remaining merger. However, 
the high gas fraction galaxies show a weak response to their initial 
pericentric passage, with strong Alog(SFR) only appearing in the 
final stages of the interaction. 

We lastly plot the evolution of the gas fraction in these galax¬ 
ies, relative to their isolated counterpart. For each time step, we 
divide the gas fraction of the pair by the gas fraction of the iso¬ 
lated galaxy. By this definition, a value of 1.0 indicates that the two 
galaxies are consuming their gas at the same rate, which is the case 
before the first close passage. As the pair galaxies consume their 
gas at a higher rate than they would in isolation in order to sustain a 
higher SFR, the ratio of the gas fractions declines. We note that for 
low gas fractions, the ratio drops more strongly than for high gas 
fractions; this is consistent with the third panel, which indicated 
that the low gas fraction galaxies had undergone a more significant 
change to their SFR, relative to their isolated counterpart, than the 
high gas fraction galaxies had. For moderate gas fractions, the gas 
consumption is expected to be rather small; [Fertig et al.| ( |2015^ find 
no evidence for large changes in gas fraction in their pairs sample, 
consistent with the results of prior simulations. 


^ Two out of the 17 data points’ removal results in a probability of null 

correlation that cannot be excluded at > 2cr. The least constraining null ® The second and third panels are both sampled at fixed time steps of 20 

probability is 0.0621, and the most 0.0030. million years. 
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Figure 3. Results from our suite of 9 galaxy merger simulations. In all panels, colours indicate the initial gas fraction, labeled in the top panel. In the top panel, 
the three dimensional separation between the two galaxies is plotted as a function of time. The galaxies all begin on the same orbit, so their tracks overlap until 
the final stages of the merger. We truncate the data plotted here once the galaxies have a final sepai'ation of 0 kpc. In the second panel, we plot the SFR for 
each galaxy as a function of time. The third panel shows the Alog(SFR) from the simulations, normalised to the SFR of the isolated galaxy of the same initial 
gas fraction, as a function of time. In the lowest panel, we plot the ratio of the gas fraction of the interacting pair, relative to the gas fraction of the isolated 
galaxy at that time step. The horizontal dashed lines indicate where the interacting pairs and the isolated control are the same. 
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Figure [^indicates that the highest gas fraction galaxies show 
higher SFRs prior to an interaction than the lowest gas fraction 
galaxies. In spite of the strong star formation prior to an interaction, 
high gas fraction galaxies are not strongly enhanced after a close 
passage. By contrast, the low gas fraction galaxies show strong en¬ 
hancements in their SFR, reaching enhancements of nearly a factor 
of 10 over their isolated counterparts shortly after a close encounter. 
The MgasIM, = 0.04 and Mg^JM, = 0.09 runs are dramatically el¬ 
evated at all times after first pericentre. As the galaxies approach 
coalescence, the highest gas fraction galaxies now show a very 
strong response to coalescence, reaching SFR enhancements of a 
factor of 30 above their isolated controls (Alog(SFR) = 1.5). The 
lowest gas fraction galaxies show a much weaker response to coa¬ 
lescence. Although all gas fractions show enhanced star formation 
after first pericentre, the relative strength of the low gas fraction 
starburst after first pericentre seems to indicate that high gas frac¬ 
tions are less efficient at forming new stars between first and second 
pericentre, but when the torques are maximal during coalescence, 
the larger gas volume may be able to create a larger burst of star 
formation. 

To compare more closely with the observations, we take two 
further steps with the simulations. First, we convert the real sepa¬ 
rations into projected separations, in order to mimic the effects of 
projection present in the observations, following the methodology 
of |Scudder et al.| j2012b| l and |Patton et al.|p013^ . For each point 
in the simulations (corresponding to individual time steps), a 3 di¬ 
mensional random viewing angle is applied to the separation data. 
The 3-D angle is defined as |cos((^)|, where (j) = sirr\R) and R is 
a different random value between 0 and 1 for each data point. The 
cos(0) term accounts for a 2 dimensional spin of the viewing angle, 
and the arcsin within the (j) term accounts for the 3 dimensional 
distribution of angles across a sphere. This is repeated for each 
time step, which produces a new random viewing angle for each 
of the data points in the simulation. This conversion from real into 
projected separations preferentially scatters data points to smaller 
separations. The projected separations are then plotted against their 
Alog(SFR) values in Figure]^ 

The galaxy pairs in our VLA sample are explicitly limited 
to r,, < 30 kpc. However, the angular separation required for 
the VLA to spatially resolve the sample imposed an indirect lower 
limit on r^, such that all galaxies in the sample fall in a relatively 
narrow range of Vp: 15 < Vp < 30 kpc. To mimic this con¬ 
straint in the simulations, we select out the subsample of the simu¬ 
lation data which falls within projected separations of between 15 
and 30 kpc, indicated in Figure]^ as the grayscale region. As 
the observations are also limited to instantaneous gas fractions, we 
plot in Figurej^the instantaneous gas fraction of the galaxy versus 
Alog(SFR), for those time steps which placed the data points within 
15 < Tp < 30 kpc. 

Figure|^illustrates Alog(SFR) as a function of log(Mjoj/M,). 
As this region of rp contains both galaxies which have not yet had 
a close pass, those just after their close pass, and galaxies returning 
towards coalescence, there is a spread in Alog(SFR) values for any 
given value of Mg^sIM,. However, it is clear that those galaxies in 
the lowest MgaJM, regime are also those galaxies which permit 
the strongest Alog(SFR) enhancements. In comparing Figure!^ to 
the results found based upon the observational data in Figureplthe 
correlations of the simulations and the observations appear to be in 
contrary directions. We discuss this apparent contradiction further 
in the Discussion section. 



Figure 4. Alog(SFR) as a function of rp. As in Figurecolours indicate 
the initial gas fraction. Each data point from the third panel of Figure[^has 
had a random viewing angle applied to its separation to imitate the effects 
of observational projection effects, which tends to make galaxies appear 
closer in the sky than their true physical separation. The grey shaded region 
indicates the region spanned by the observational sample. 



Figure 5. log(Mjaj/M,) as a function of Alog(SFR), for points which fall 
between 15 and 30 /tjg* kpc after projection effects. Colours indicate the ini¬ 
tial gas fraction of the simulation, shown in the legend. Galaxies with the 
highest instantaneous gas fractions are not typically found to have strong 
Alog(SFR) enhancements within this range of rp. However, the galaxies 
with low instantaneous gas fractions show a wide range of Alog(SFR), 
which include enhancements up to ~ 0.9. 
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6 DISCUSSION 

The current observational results suggest a positive correlation be¬ 
tween the HI content of a galaxy and the strength of triggered star 
formation within that galaxy. Spearman rank correlation tests in¬ 
dicate that the data exclude the null hypothesis at ~ 2.5cr signif¬ 
icance, and a jackknife resampling of the data indicates that the 
correlation is not extremely sensitive to individual galaxies in the 
sample. In contrast, simulations of binary mergers with varying gas 
fractions indicate that prior to coalescence, we would expect an 
anti-correlation between gas fraction and SFR enhancement, rela¬ 
tive to an identical galaxy in isolation. We now wish to investigate 
what may lie at the root of this difference. 


6.1 Reconciling the different physical pictures 

There is a critical difference between the simulations and the cur¬ 
rent observational data; for the simulations, the control is an exact 
match in all parameters to the interacting galaxy prior to its interac¬ 
tion. The observational data is very tightly matched in stellar mass, 
redshift, and local environment, but, importantly, the gas fractions 
of the SDSS control galaxies are unconstrained. 

We test the impact of leaving the gas fraction unconstrained 
in the observational controls using the simulation data. We choose 
a single simulated galaxy (in isolation) to operate as the ‘average’ 
control galaxy for all of the interacting galaxies. We then recal¬ 
culate the Alog(SFR) using this fixed control galaxy for all simu¬ 
lated encounters. This means that every interacting galaxy’s SFR 
is being compared to the same isolated galaxy, which has the same 
stellar mass but a different gas fraction. For this recalculation of 
Alog(SFR), we select the Mg^sIM, = 0.32 run (an intermediate 
gas fraction) as our ‘average’ control. In all other aspects, the cal¬ 
culation proceeds as before. We regenerate Figure ^with the re¬ 
calculated Alog(SFR) values in Figure]^ Figure |6 shows the in¬ 
stantaneous gas fraction as a function of Alog(SFR) for galaxies 
within 15 <r,,< 30 kpc, now using the fixed control galaxy 
MgaJM, = 0.32 for the calculation of Alog(SFR). 

A comparison of Figure and Figure demonstrates that 
the lack of a gas fraction match in the control sample drasti¬ 
cally changes the apparent relationship between gas fraction and 
Alog(SFR) from an anti-correlation to a correlation. By using a 
galaxy with an average gas fraction as the simulation’s control, the 
trend between gas fraction and Alog(SFR) in the simulations is now 
in broad agreement with the observations. As demonstrated in the 
second panel of Figure]^ there is a strong pre-interaction correla¬ 
tion between gas fraction and SFR, which can explain this change. 
Although galaxies with low initial gas fractions exhibit the highest 
interaction-triggered enhancements, they begin an interaction with 
low absolute SFRs. These initially low absolute SFRs translate to 
strongly negative Alog(SFR)s when compared with an average gas 
fraction control. For these low gas fraction systems, even the strong 
increase in their SFRs during the interaction is still not enough of 
a boost to increase their SFRs above the non-interacting SFR of an 
average gas fraction control. Conversely, high gas fraction galaxies 
have a pre-interaction SFR that is already enhanced relative to the 
average gas fraction control. Although their interaction-induced en¬ 
hancement is modest, they have such strong SFR at the outset that 
they still exhibit the highest Alog(SFR)s. For galaxies without a 
control sample explicitly matched in gas fraction, the correlation 
between the gas fraction of a galaxy and its SFR, independent of 
the interaction environment, is clearly much stronger than the cor- 



Total Alog(SFR) 

Figure 6. A counterpart to Figure|^ with Alog(SFR) values calculated rel¬ 
ative to a fixed gas fraction isolated galaxy {Mgas/M,= 0.32). There now 
appears to be a strong correlation between the initial gas fraction and the 
level of SFR enhancement. 


relation between the gas fraction and the triggered SFR due to an 
interaction. 

In the ideal case, this test would lead us to construct a more 
suitable control sample for our observational sample in order to re¬ 
move this source of bias. However, we would require a large sam¬ 
ple with appropriately deep H i observations in order to construct 
a control sample matched in gas fraction such a sample does not 
currently exist. ALFALFA.40 ( [Haynes et al.|[20TT| > is the largest 
sample of galaxies with cross-matches to the SDSS, which would 
normally produce the best number statistics for the construction of 
a control sample. However, our VLA sample is significantly deeper 
than the ALFALFA.40 sample at a given stellar mass. This differ¬ 
ence in survey depth means that we would only be able to find gas 
fraction matches to the highest gas fraction galaxies at fixed stellar 
mass, and these are the galaxies that the simulations indicate are 
least likely to show a strong enhancement. Without the ability to 
match our entire VLA sample, the ALFALFA sample of galaxies is 
not useful as a source of control galaxies. 

We note that the simulation suite presented here provides a 
wide range of gas fractions at fixed stellar mass, and therefore 
should serve as an upper limit to the strength of correlation driven 
by poor matching in gas fraction. Even though log(MgaJM,) spans 
more than an order of magnitude in the simulations, a similar range 
of HI gas fractions at fixed stellar mass has been observed in large 
Hi surveys (e.g., jCatinella et al.|2012[|Huang et al.||2012[ , indi¬ 
cating that simple stellar mass matching is unlikely to indirectly 
control for the gas fraction of the galaxy, and this source of bias is 
almost certainly present in the observational sample. 
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6.2 The Impact of Molecular Gas 


As described in Section [4T| our adopted definition of gas fraction 
in this work does not include any contribution from the molecular 
gas component of a galaxy. However, the impact of the exclusion 
of molecular gas is not well understood. |Saintonge et al.| |20TT| 
finds that for non-interacting galaxies, the molecular gas mass of a 
galaxy is about 30% of the neutral gas mass on average, with a great 
deal of scatter; the two quantities are, at best, only weakly corre¬ 
lated. The behaviour of the molecular gas in interacting galaxies is 
still very unclear. It has been suggested that in galaxy mergers, the 
relationship between SFR and the molecular gas content is system¬ 
atically changed such that interacting galaxies have significantly 
faster depletion times; i.e., they turn more gas into stars per unit 
time, when compared to their isolated counterparts (|Daddi et al. 


|20^ Genzel et al.|2010{ [Saintonge et al.|2012| >, with Stark et al. 


12013 I finding that recently enhanced SFR (such as is present in in¬ 
teractions) correlates with a larger H 2 /H 1 content. However, |Zasov| 
|& Sulen tT^ jl994^ finds that interacting galaxies have only a small 
component of molecular gas, and |Casasola et al!] ( |2004^ found that 
galaxies in pairs have normal H i content but higher than average 
molecular gas (resulting in a higher than average total gas mass for 
interacting galaxies). Until a direct test of the dependence of SFR 
on the H 2 content of a sample of galaxy pairs is undertaken, it will 
remain unclear how the molecular gas will correlate with the SFR 
enhancement of a galaxy. 

The measurements from the simulations presented here (to 
which we are comparing the observations) are akin to the the to¬ 
tal gas mass of the galaxy, with some fraction of the total gas mass 
divided between Hi and H 2 . While it would be interesting to ex¬ 
amine the impact of molecular gas in simulated interactions, this is 
beyond the scope of the current work. 


7 CONCLUSIONS 

In this work, we have attempted to determine if the gas content of 
a galaxy plays a dominant role in influencing the strength of an 
interaction-triggered starburst. To tbis end, we present new obser¬ 
vational data from the VLA, which has allowed us to estimate the 
Hi content of a sample of 34 galaxies in 17 galaxy pairs. Of the 
34 galaxies, 17 are detected with S/N > 3.0. The detected Hi mass 
contained within these galaxies was converted into a gas fraction, 
using the stellar masses of |Mendel et al.|p014^ . 

We compare the gas fraction of these galaxies with their SFR 
enhancements (Alog(SFR)), which are determined relative to a con¬ 
trol sample which is tightly matched in stellar mass, redshift, and 
local density. We find that the H i gas fraction of the VLA galaxies 
is correlated at ~ 2.5cr with the total (i.e., across the entire galaxy) 
star formation enhancement of that galaxy. We then compare the 
observational results with a suite of simulations in order to help 
develop an interpretive framework for the observations. 

The simulation suite consists of a set of 9 simulated galaxy 
mergers which vary only in gas fraction. The simulations suggest 
that the lowest gas fraction galaxies should show the highest en¬ 
hancements to their SFR after their first encounter with their com¬ 
panion. High gas fraction galaxies appear the least affected by their 
close pass, relative to an identical, non-interacting galaxy, implying 
an anti-correlation between Alog(SFR) and gas fraction. 

The simulations also indicate that high gas fraction galaxies 
have systematically higher pre-interaction SFRs when compared to 
low gas fraction galaxies. We determine that if the gas fractions of 


the control galaxies are not part of the matching procedure (as is the 
case for our observational sample), an apparent correlation between 
Alog(SFR) and gas fraction can be produced, as the correlation be¬ 
tween higher base SFR and higher gas fraction is much stronger 
than the anti-correlation of Alog(SFR) and gas fraction due to the 
interaction itself. If the simulations are handled in a manner consis¬ 
tent with the observations, we find the same qualitative correlation 
between gas fraction and SFR enhancement. 
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Figure Al. The log of the gas fraction vs. Alog(SFR). Points are colour- 
coded and labeled with a greek letter according to pair. Large circles indi¬ 
cate S/N > 3.0, with smaller triangles indicating the 3(T upper limits for the 
maximum amount of flux that could exist in the spectrum. 
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APPENDIX A: ADDITIONAL FIGURES 

In the Appendix, we present figures displaying the full data set. 
Figure [AT| shows the full sample of 34 galaxies in our sample, in¬ 
cluding all non-detections. Figures [A^ - | A17 | show the SDSS image 
of the galaxy pair, along with the associated H i spectrum for each 
of the galaxies in the sample. All images are at the same angular 
scale for consistency, and are 250 arcsec to a side. 
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Figure A2. Left panel: SDSS thumbnail of galaxy pair 587727178473930875 (right) & 587727178473930886 (left). Right panel: Flux density vs. optical 
velocity for galaxy pair 587727178473930875 (upper) & 587727178473930886 (lower). The vertical dashed line indicates the velocity of the galaxy, horizontal 
dotted line indicates 0 flux. The blue shaded region indicates the region of the spectra used for signal to noise and gas mass calculations. The grey shaded 
region around the indicated velocity of the galaxy indicates the uncertainty of the SDSS redshift (±30 km s“^). The upper spectrum has a peak/RMS S/N of 
2.62, while the lower spectrum has a peak/RMS S/N of 2.24. 
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Figure A3. Same as in Figure [A^ Left panel: SDSS thumbnail of galaxy pair 588017702411763744 (lower left) & 588017702411763872 (upper right). Right 
panel: Flux density vs. optical velocity for galaxy pair 588017702411763744 (upper) & 588017702411763872 (lower). The upper spectrum has a peak/RMS 
S/N of 2.67, while the lower spectrum has a peak/RMS S/N of 2.79. 
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HI gas in galaxy pairs 17 




Figure A4. Same as in Figure fA^ Left panel: SDSS thumbnail of galaxy pair 587739303684866183 (upper left) & 587739303684866173 (lower right). Right 
panel: Flux density vs. optical velocity for galaxy pair 587739303684866183 (upper) & 587739303684866173 (lower). The upper spectrum has a peak/RMS 
S/N of 2.35, while the lower spectrum has a peak/RMS S/N of 3.57. 
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Figure AS. Same as in Figure [A^ Left panel: SDSS thumbnail of galaxy pair 588017605758025795 (lower) & 588017605758025732 (upper). Right panel: 
Flux density vs. optical velocity for galaxy pair 588017605758025795 (upper) & 588017605758025732 (lower). The upper spectrum has a peak/RMS S/N of 
1.93, while the lower spectrum has a peak/RMS S/N of 3.62. 
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18 Scudder et al. 



Figure A6. Same as in Figure [A^ Left panel: SDSS thumbnail of galaxy pair 587742901789589569 (right) & 587742901789589575 (left). Right panel: Flux 
density vs. optical velocity for galaxy pair 587742901789589569 (upper) & 587742901789589575 (lower). The upper spectrum has a peak/RMS S/N of 3.20, 
while the lower spectrum has a peak/RMS S/N of 2.70. 
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Figure A7. Same as in Figure [A^ Left panel: SDSS thumbnail of galaxy pair 588023670245949622 (upper) & 588023670245949625 (lower). Right panel: 
Flux density vs. optical velocity for galaxy pair 588023670245949622 (upper) & 588023670245949625 (lower). The upper spectrum has a peak/RMS S/N of 
1.61, while the lower spectrum has a peak/RMS S/N of 2.73. 
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HI gas in galaxy pairs 19 





Figure A8. Same as in Figure fA^ Left panel: SDSS thumbnail of galaxy pair 587733605328093368 (upper right) & 587733605328093256 (lower left). Right 
panel: Flux density vs. optical velocity for galaxy pair 587733605328093368 (upper) & 587733605328093256 (lower). The upper spectrum has a peak/RMS 
S/N of 2.26, while the lower spectrum has a peak/RMS S/N of 4.04. 



Figure A9. Same as in FigurejA^ Left panel: SDSS thumbnail of galaxy pair 587727179536859247 (upper left) & 587727179536859227 (lower right). Right 
panel: Flux density vs. optical velocity for galaxy pair 587727179536859247 (upper) & 587727179536859227 (lower). The upper spectrum has a peak/RMS 
S/N of 11.40, while the lower spectrum has a peak/RMS S/N of 2.02. 
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20 Scudder et al. 



Figure AlO. Same as in Figure [A^ Left panel: SDSS thumbnail of galaxy pair 587729158970867777 (right) & 587729158970867792 (left). Right panel: 
Flux density vs. optical velocity for galaxy pair 587729158970867777 (upper) & 587729158970867792 (lower). The upper spectrum has a peak/RMS S/N of 
4.09, while the lower spectrum has a peak/RMS S/N of 3.71. 





Figure All. Same as in Figure [A^ Left panel: SDSS thumbnail of galaxy pair 587739609695453284 (upper) & 587739609695453281 (lower). Right panel: 
Flux density vs. optical velocity for galaxy pair 587739609695453284 (upper) & 587739609695453281 (lower). The upper spectrum has a peak/RMS S/N of 
4.08, while the lower spectrum has a peak/RMS S/N of 3.69. 
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Hi gas in galaxy pairs 21 
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Figure A12. Same as in Figure [A2| Left panel: SDSS thumbnail of galaxy pair 588848899908370674 (lower) & 588848899908370505 (upper). Right panel: 
Flux density vs. optical velocity for galaxy pair 588848899908370674 (upper) & 588848899908370505 (lower). The upper spectrum has a peak/RMS S/N of 
4.54, while the lower spectrum has a peak/RMS S/N of 3.26. 



Figure A13. Same as in Figure [A^ Left panel: SDSS thumbnail of galaxy pair 587726033308680234 (lower) & 587726033308680320 (upper). Right panel: 
Flux density vs. optical velocity for galaxy pair 587726033308680234 (upper) & 587726033308680320 (lower). The upper spectrum has a peak/RMS S/N of 
1.97, while the lower spectrum has a peak/RMS S/N of 2.00. 
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Figure A14. Same as in Figure [A^ Left panel: SDSS thumbnail of galaxy pair 587741489815027774 (right) & 587741489815028146 (left). Right panel: 
Flux density vs. optical velocity for galaxy pair 587741489815027774 (upper) & 587741489815028146 (lower). The upper spectrum has a peak/RMS S/N of 
2.70, while the lower spectrum has a peak/RMS S/N of 2.23. 





Figure A15. Same as in Figure [A2| Left panel: SDSS thumbnail of galaxy pair 587744873717563559 (right) & 587744873717563471 (left). Right panel: 
Flux density vs. optical velocity for galaxy pair 587744873717563559 (upper) & 587744873717563471 (lower). The upper spectrum has a peak/RMS S/N of 
5.35, while the lower spectrum has a peak/RMS S/N of 3.51. 
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Figure A16. Same as in Figure fA^ Left panel: SDSS thumbnail of galaxy pair 587729160043757697 (lower right) & 587729160043757707 (upper left). Right 
panel: Flux density vs. optical velocity for galaxy pair 587729160043757697 (upper) & 587729160043757707 (lower). The upper spectrum has a peak/RMS 
S/N of 3.76, while the lower spectrum has a peak/RMS S/N of 3.16. 
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Figure A17. Same as in FigurejA^ Left panel: SDSS thumbnail of galaxy pair 587726033341776175 (upper right) & 587726033341776191 (lower left). Right 
panel: Flux density vs. optical velocity for galaxy pair 587726033341776175 (upper) & 587726033341776191 (lower). The upper spectrum has a peak/RMS 
S/N of 2.85, while the lower spectrum has a peak/RMS S/N of 3.42. 
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